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Abstract 

The NA60 experiment has studied J/ip production in p-A collisions at 158 
and 400 GeV, at the CERN SPS. Nuclear effects on the J/ip yield have 
been estimated from the A-dependence of the production cross section ratios 
a f/il>l a fU (A=A1, Cu, In, W, Pb, U). We observe a significant nuclear sup- 
pression of the J/ip yield per nucleon-nucleon collision, with a larger effect 
at lower incident energy, and we compare this result with previous obser- 
vations by other fixed-target experiments. An attempt to disentangle the 
different contributions to the observed suppression has been carried out by 
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studying the dependence of nuclear effects on x 2 , the fraction of the nucleon 
momentum carried by the interacting parton in the target nucleus. 
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1. Introduction 

The study of charmonium production in hadronic collisions is an interest- 
ing test of our understanding of the physics of strong interactions. While the 
production of the cc pair can be addressed in a perturbative-QCD approach, 
the subsequent binding of the pair is an essentially non-perturbative process, 
involving soft partons and occurring on a rather long timescale (> 1 fm/c). 
Various theoretical approaches have been proposed (see [l[ for a recent re- 
view), but a satisfactory description of J/0 production in p-p collisions is 
still missing jij. In p-A interactions, the heavy-quark pair is created in the 
nuclear medium, and the study of its evolution towards a bound state can 
add significant constraints to the models. For example, the strength of the 
interaction between the evolving cc pair and the target nucleons, that can 
lead to a break-up of the pair and consequently to a suppression of the J/ip 
yield, may depend on its quantum state at the production level (color-octet 
or color-singlet), and on the kinematic variables of the pair jsl, 0). In addi- 
tion to final-state effects, also initial-state effects may influence the observed 
J /if) yield in p-A. In particular, parton shadowing in the target nucleus [i[ 
may suppress (or enhance, in the case of anti-shadowing) the probability of 
producing a J/0, while the energy loss of the incident parton [f| in the nu- 
clear medium, prior to cc production, may significantly alter the J/ip cross 
section and kinematic distributions. Furthermore, effects as final-state en- 
ergy loss and the presence of an intrinsic charm component in the proton 
may also play a significant role 0]. Clearly, the correct understanding and 
the disentangling of the various nucleus-related effects on J/0 production is 
a non-trivial task, which poses significant challenges to theory, but at the 
same time offers important insights on the J/ip production and interaction 
mechanisms. Finally, a suppression of the J/ip has been proposed a long 
time ago as a signature of the formation, in ultrarelativistic nucleus-nucleus 
collisions, of a state where quarks and gluons are deconfined (Quark-Gluon 
Plasma) Results from p-A collisions, taken in the same kinematic con- 
ditions of A-A, and properly extrapolated to nucleus-nucleus collisions, are 
therefore necessary to calibrate the contribution of the various cold nuclear 
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matter effects to the overall observed suppression 0,0. 

Having to deal with a rather complicated interplay of various physical 
processes, the availability of accurate sets of data, spanning large intervals in 
the incident proton energy, and covering large xf and px regions, is essential 
for a thorough understanding of the involved mechanisms. At fixed target 
energies, high-statistics J/ip samples have been collected in recent years by 



the DESY experiment HERA-B [ill], at 920 GeV incident energy, by E866 [12 



at FNAL at 800 GeV, and by the CERN-SPS experiment NA50 at 400 and 



450 GeV 13 



In this Letter, we present an extension of these measurements towards 



lower energies, carried out by the NA60 experiment [14j. J/0 production 
has been studied in p-A collisions at 158 GeV, the same energy used in 
A-A collisions at the CERN SPS. Data have also been taken with a 400 
GeV beam, in order to provide a result that can be compared with previous 



sets of p-A data taken by the NA50 experiment [13J. We have performed a 
systematic study of nuclear effects by analyzing the A— dependence of the 
J/i/j production cross section on seven different target nuclei (Be, Al, Cu, In, 
W, Pb and U). The results are relative to a region close to midrapidity, and 
are presented differentially in zp and £2 and compared with previous results 
from the higher- -y/i measurements mentioned above. The influence of parton 
shadowing in the nuclear targets on the results is also discussed, based on 
various recent parametrizations of this effect (5), [jjl, 16, 17]. A comparison 
between our results at the two energies as a function of X2, the fraction of the 
nucleon momentum carried by the parton in the nuclear target that produces 
the J/0, is carried out. By doing so, the shadowing effects can be factorized 
and the importance of the other nuclear effects can be investigated. 

2. Data analysis 

The NA60 experiment has measured muon pair production in p-A and 
A-A collisions at the CERN SPS. Its experimental apparatus was based on 
a muon spectrometer (MS), positioned downstream of a hadron absorber 
with a total thickness of 12 nuclear interaction lengths (A/). The MS is 
coupled to a vertex spectrometer (VT) based on Si pixel detectors. The 
experiment triggered on muon pairs detected in the MS, which were then 
matched, during reconstruction, to the corresponding tracks in the VT. For 
details on the detector set-up and matching between the MS and VT we refer 

to 05 
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The J/ijj mesons have been identified through their decay to a muon pair. 
In the p-A data taking, the SPS proton beam, with an average intensity of 
5-10 8 s _1 , was hitting a target system composed of nine sub-targets with 
thicknesses between 0.005 and 0.012 A/, with relative spacing between tar- 
gets from 0.8 to 1 cm. To limit the possible influence of position-dependent 
systematic effects on the measurement of the nuclear dependence of the J/i/j 
yield, the subtargets were placed in a mixed-A order, namely Al, U, W, Cu, 
In, Be, Be, Be and Pb. 

The analysis described in this Letter has been carried out on a data sample 
consisting of 2.8-10 6 events at 158 GeV, and 1.5-10 6 events at 400 GeV, 
containing a dimuon reconstructed in the MS. The mass region m MM >2.85 
GeV/c 2 , dominated by J/ip decays, contains 3.2-10 4 and 2.1-10 4 events at 
158 and 400 GeV, respectively. For about 50% of these events, each of the 
two MS tracks can be matched in direction and momentum with a track in 
the VT. Such a value of the dimuon matching efficiency is due the smaller 
solid angle coverage of the VT with respect to the MS and to the inefficiency 
of the pixel sensors as will be discussed below. A Monte-Carlo simulation 
has shown that the contamination from events where at least one of the MS 
tracks is matched to a wrong VT track is negligible in the J/ip mass region. 

For events containing a pair of hard muons from the decay of a J/ip (the 
average muon momentum is 27 GeV/c), the point of closest approach of 
the two muons gives an accurate estimate of their production point, with a 
resolution of ~ 650 fim. In Fig. [T]we show the distribution of the longitudinal 
coordinates of the dimuon production vertex, for m MA1 >2.85 GeV/c 2 . The 
peaks corresponding to the position of the various production targets are 
clearly visible and well separated, showing that the nuclear target where the 
J/ip has been produced is unambiguously determined. 

The MS detects J/ip produced in the rapidity range 3< yi a b <4. However, 
the acceptance of the VT for matched tracks from the J/ip decay is target- 
dependent, the more downstream targets covering smaller rapidities. We 
have therefore selected events produced in the range 3.2< yi^ <3.7, which is 
covered for all the targets. The cut -0.5< cos 6^3 <0.5, on the polar angle of 
the decay muons in the Collins-Soper reference frame has also been applied, 
in order to remove events at the edge of the MS's acceptance. The transverse 
momentum coverage for the J/ijj measurement extends down to zero pt- The 
acceptance varies by about 30% in the pt range accessible with the collected 
statistics (px ^ 3 GeV/c). 

The J/ip signal has been extracted for each target nucleus at the two 
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Figure 1: The distribution of the longitudinal coordinate of the point of closest approach 
of the opposite sign muon pair, for m w >2.85 GeV/c 2 . 



energies by fitting the opposite sign mass spectra, in the region >1.5 
GeV/c 2 , with a superposition of the mass shapes of the expected dimuon 
sources. These include the Drell-Yan process (DY), the semi-leptonic decays 
of correlated D- meson pairs (DD), the J/ip and the ip' resonances. The 
expected mass distributions for the DY and DD contributions have been 



calculated with PYTHIA |l8j, using the GRV94LO [19| parton distribution 



functions. The J/^ events have been generated, for the 400 GeV data, using 



the y and distributions measured with good accuracy by NA50 20]. At 



158 GeV, the differential distributions have been tuned directly on the data, 
using an iterative procedure. The events have been tracked through the set- 
up and then reconstructed with the same algorithm used for real data. The 
contribution from the combinatorial pair background due to it and K decays 
(completely negligible in the J/ijj mass region) has been estimated through an 
event mixing technique 14j. In Fig. [2] we show, as an example, the opposite- 
sign dimuon invariant mass distribution, relative to the kinematical domain 
specified above, for p-In interactions at 158 GeV. The quality of all the fits 
to the invariant mass spectra is satisfactory (with x 2 /ndf ranging from 0.7 to 
1.5). The mass resolution at the J/ip peak is ~70 MeV/c 2 , and the number 
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of J/ip events ranges from ~800 to ~2000, depending on the target. The 
number of continuum events under the J/ip peak is very small (<4%), and 
is dominated by the DY contribution. 



X 2 / ndf 82.98 / 65 




m m (GeV/c 2 ) 

Figure 2: Fit to the /x + ' [i~ invariant mass spectrum for p-In collisions at 158 GeV. The 
dashed line represents the Drell-Yan process, the dotted lines the charmonium resonances, 
the dashed-dotted line the DD contribution, the thin continuous line the combinatorial 
background. The thick continuous line is the sum of all the contributions. 

Nuclear effects have been parametrized by fitting the A— dependence of 
the production cross section with the simple power law (Xj^ = • A a , and 
then studying the evolution of a as a function of various kinematic variables. 
Alternatively, nuclear effects have been quantified by fitting the data in the 
framework of the Glauber model, having as input parameters the density 
distributions for the various nuclei. The model gives as output the so-called 
J/ip absorption cross section crfj^. Clearly, both a and crfj s ^ represent effective 
quantities, including the contribution of the various sources of nuclear effects 
detailed in the Introduction. 

The nuclear effects on J/?/) production have been evaluated starting from 
the cross section measured for each target, normalized to the cross section 
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for the lightest one (Be): 

<t 3 1+/a _ N^/A I N%* I ' A Be 

4t I A Be = NT ■ NT S -A A -e A / iVg • N^f ■ A Be ■ e Be W 

where, for the target with mass number A, is the number of J/0 events, 
N A nc is the number of incident protons, iV j4 arg is the number of target nuclei 
per unit surface, A a is the J/-0 acceptance and is the detection efficiency. 
When building these ratios, the results obtained with the three Be targets 
have been averaged. The use of relative cross sections brings several advan- 
tages. In particular, the number of incident protons cancels out (a small 
attenuation factor, due to the inelastic scattering of beam particles, which 
is 6% for the most downstream target, has been corrected for) since all the 
targets were simultaneously exposed to the beam, which had a transverse 
dimension much smaller than that of the targets. Furthermore, the fraction 
of the detection efficiency related to the MS also cancels out. In fact, this 
detector cannot distinguish the target where the dimuon has been produced. 
This is due to the presence of the thick hadron absorber, to the large dis- 
tance between the target system and the MS tracking chambers (6 to 16 m), 
and to the closely spaced targets. Therefore, the muon detection efficiency 
is independent of the production target. 

Contrary to the MS case, the angular acceptance covered by the VT is 
slightly target dependent. Therefore, non- uniformities in the efficiency of a 
certain pixel plane lead to a rapidity-dependent inefficiency which is different 
from target to target. The efficiency of each of the 15 VT planes has been 
estimated using a modified track reconstruction algorithm that excludes the 
plane under study. For each reconstructed track we then check the presence of 
a hit in that plane in a fiducial region around the intersection of the track with 
the plane. Efficiency values are calculated on a run-per-run basis (<2 hours), 
for sub-regions of the pixel plane down to a size of ~0.2 mm 2 , if the track 
statistics in the sub- region under study is large enough (more than 40 tracks). 
Otherwise, contiguous regions are grouped in order to reach a statistically 
significant track sample. The distribution of the efficiency values for the 
various regions of the VT is peaked at ~90%, with 82% of the detector area 
having an efficiency larger than 60%. Finally, the J/ip acceptances, estimated 
by means of the Monte-Carlo simulation, range between 14.9% and 23.0% 
in the kinematic range 3.2< y [ab <3.7, |cos#cs| <0.5. The quoted rapidity 
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Table 1: The relative values, target by target, of the integrated luminosities, and of the 
products Aa ■ £A- We have chosen as reference value the most upstream Be target. The 
Al target was not in place during the data taking at 400 GeV. 



range corresponds to the center-of-mass rapidity windows 0.28< y <0.78 at 
158 GeV and -0.17 < y <0.33 at 400 GeV. 

In Table 1 we summarize, target by target and for the two energies under 
study, the relative values of the J/^ acceptance times detection efficien- 
cies. For completeness, we also report the relative integrated luminosities 
C A /C Be = (NT ■ NT g )/(N™ • N^). 



3. Results 

In Fig. [3] we present the cross section ratios (cr/^ / Aj) / (cr J J^ / A Be ) , at 
158 and 400 GeV, where A t is the nuclear mass number of target i. The 
results are shown as a function of L, the mean thickness of nuclear mat- 
ter crossed by the J/ip in its way through the target nucleus. The L val- 
ues have been computed with the Glauber model, using realistic density 



distributions for the various nuclei [2l|. The quoted systematic uncertain 



ties include contributions, quadratically combined, coming from the uncer- 
tainty on i) the measurement of the target thicknesses (<1.5%), ii) the 
J/ip acceptance, due to the choice of the rapidity distribution adopted in 
the Monte-Carlo calculation (<1.5%), hi) the efficiency calculation (<3.0%). 
This last, and most important, contribution has been obtained varying by 
±10% the estimated efficiency values of the VT pixel detectors. We only 
quote, for each incident energy, the fraction of the systematic uncertainty 
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which is not common to all the points, the only one relevant when plot- 
ting relative cross sections. Fig. [3] shows, for both datasets, a suppression 
of the J/0 yield when moving from light to heavy targets, and, in particu- 
lar, a larger suppression for the 158 GeV data sample. Using the Glauber 
model, we have estimated <7j^(158GeV) = 7.6 ± 0.7(stat) ± 0.6(syst) mb 
and 0-^(400 GeV) = 4.3 ± 0.8(stat) ± 0.6(syst) mb. With the power-law 
parameterization we get a J/v ,(158 GeV) = 0.882 ± 0.009(stat) ± 0.008(syst) 
and aj/^(400GeV) = 0.927 ± 0.013(stat) ± 0.009(syst). We note that the 
400 GeV values are in good agreement with results from NA50 at the same 
energy (cij^ = 4.6 ± 0.6 mb [13J), in a similar rapidity range. We also stress 
that the reported values have not been corrected for shadowing effects and 
therefore represent a global effective estimate of the nuclear effects on the 
J/0 yield. 
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Figure 3: Cross sections for J/-0 production in p-A collisions, normalized to the p-Be 3/tp 
cross section. The squares represent the 158 GeV results, the circles refer to 400 GeV. 
The lines are the fit results. 

To get further insight on the dependence of nuclear effects on the J/?/' 
kinematic variables we compare in Fig. HJ as a function of xp, the a values 
obtained in this analysis with those from previous experiments in the fixed- 
target energy range. For this purpose, the NA60 sample has been subdivided 
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Energy (GeV) 


Xp 


a 


158 


05 — 15 


911 ± 019 ± 036 


158 


0.15^0.20 


0.868 ±0.016 ±0.009 


158 


0.20^0.25 


0.887 ±0.017 ±0.008 


158 


0.25^0.30 


0.867 ±0.019 ±0.010 


158 


0.30 -r- 0.40 


0.856 ±0.025 ±0.027 


400 


-0.075^-0.025 


0.924 ±0.026 ±0.021 


400 


-0.025 ^0.025 


0.920 ±0.019 ±0.012 


400 


0.025 ^0.075 


0.924 ±0.020 ±0.013 


400 


0.075 ^0.125 


0.906 ±0.028 ±0.025 



Tabic 2: a as a function of Xp for the 158 and 400 GeV data sample. The first quoted 
uncertainty is statistical, the second is the systematic one. 



into 5 xp bins at 158 GeV (4 at 400 GeV), covering the region 0.05 < xp < 
0.40 (—0.075 < xp < 0.125). In Table 2 we present the corresponding a 
values. 

Two main features emerge from this comparison. First, when going from 
negative towards positive xp, a steadily decreases. This effect was already 
known from the data of HERA-B and E866, taken at rather similar incident 
proton energies (920, 800 GeV), and a similar effect might be present in our 
results at 158 GeV, even if the size of the errors is not negligible. Second, 
at a constant xf, the a values are lower when the incident proton energy is 
smaller, as can be seen when comparing the HERA-B/E866 results with our 
results at 400 and 158 GeV. On the other hand, it is also worth mentioning 
that from NA3 results on J/i/j production at 200 GeV (22j one extracts a 
values which are in partial contradiction with these observations, being simi- 
lar to those obtained with the higher energy data samples (HERA-B/E866). 
We also note that in our calculation of the NA3 a values, performed starting 
from their measured J/ip cross section ratios between p-Pt and p-p collisions, 
the small bias 23] induced by the use of a light target in the determination 
of a has been corrected for. 

A satisfactory theoretical interpretation of the complex observed pattern 
is missing for the moment. Various works have underlined the importance of 
several effects, including final state break-up, parton shadowing, initial and 
final state energy loss, and the presence of a charm component in the nucleon 



wavefunction [7|, |24], |25|. However, the relative weight of these effects is still 
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under debate. 
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Figure 4: The xp dependence of the a parameter. Open triangles correspond to HERA-B 
results, open squares to E866, open circles to NA50 (450 GeV), closed circles to NA60 (400 
GeV), closed squares to NA60 (158 GeV), open stars to NA3. The error bars represent 
the quadratic sum of statistical and systematic uncertainties. 



The new results from NA60 at 400 and 158 GeV have been obtained with 
the same experimental apparatus and in very similar running conditions. 
Their direct comparison may therefore be a clean testing ground for models. 
Interesting information can also be obtained by a simple comparison of nu- 
clear effects as a function of various kinematic variables. As an example, we 
now consider the x<i dependence of the a parameter, rc 2 = TWr/yfs • exp(— y) 
being the fraction of the momentum of the target nucleon carried by the 
parton which produces the J/V0- This kinematic variable is particularly 



Such a relation, commonly adopted, implies that the 2— >2 kinematics of the main J/ijj 
production process gg — > J/t/jg is, effectively, that of a 2— >1 process (gg — > J/i/>) i.e. the 
final state gluon is very soft. The alternative 2— >2 approach is discussed in [26j | . 
We also note that the relation strictly holds when px -C Pl in the cm. frame. Due to the 
low \fs, this is only approximately true and induces a smearing in the estimation of x-i . 
However, in the investigated kinematic domain, the effect is similar for the two energies. 
The resulting ~ 10% shift in the X2 values does not affect the discussion in the text. 
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interesting since x 2 is the quantity that determines the amount of shad- 
owing in the target nucleus (or anti-shadowing in our kinematic range). 
Furthermore, the center of mass energy of the J/^-N system, which is a 
relevant quantity for the final state break-up of the J/0 in the cold nu- 
clear medium, is, to a very good approximation, a function of x 2 alone 

(v/Sj/^N ~ mj/tf/ y/(l+X 2 )/ X 2 ). 

Since our results at the two energies cover, with good approximation, the 
same x 2 domain, we have performed an analysis of nuclear effects in five x 2 
bins, covering the region 0.08 < x 2 < 0.14. In Fig. [5] we plot a as a function 
of x 2 for the two energies under study. We also show in the same plot the 
expected a value that would be obtained if shadowing were the only nuclear 
effect to be present. Various parameterizations of nuclear shadowing have 
been considered ji], [l5[ 16, 17], including the recent EPS09 result, where an 
estimate of the uncertainty of the calculation was also carried out. Looking 
at Fig. [5] we see that shadowing alone would lead to a values larger than 1 in 
our x 2 acceptance. This result implies that the other nuclear effects produce 
a significant suppression of the J/-0 yield. 

The other remarkable feature emerging from Fig. [5] is that a is not the 
same at a fixed x 2 , for the two values of s/s. In our analysis, the system- 
atic uncertainties related to the measurement of a are partly correlated, 
in particular those corresponding to the measurement of the J/ip detection 
efficiency Therefore, in Fig. we plot, as a function of x 2 , the quantity 
Aa = a(400GeV) — a(158GeV) which is affected by a significantly smaller 
systematic uncertainty. The results clearly indicate that Aa ^ (in the cov- 
ered x 2 range we have (Aa) = 0.056 ± O.Oll(stat) ± 0.004(syst)), and that 
Aa does not vary appreciably in the x 2 region under study. 

Having factorized the effect of shadowing, which is the same at the two 
energies at fixed x 2 , the observation of Aa > could be attributed to a 
stronger final state interaction at 158 GeV. Since, for a given x 2 , the kine- 
matics of the collision of the J/ip with the target nucleons is the same at 
158 and 400 GeV, the observed effect may be related to a change in the J/ijj 
break-up cross section at fixed y/sJJ^N. A different weight of the color-octet 
and color-singlet precursor cc states in the production process between the 
two energies could be a natural explanation. However, Non-Relativistic QCD 
(NRQCD) calculations carried out in the fixed target energy domain 271 ]. 
only show a weak y/s dependence of the color-octet and color-singlet rela- 
tive contributions. Therefore, a strong variation of the J/i/j break-up cross 
section with y/s seems unlikely. In such a case one should observe, for the 
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two incident proton energies, the same a at a given x 2 m contrast with our 
observation. In this scenario, processes different from final state break-up 
and parton shadowing must be advocated in order to explain the observed 
nuclear effects on J/ip production. In particular, effects such as the initial 
state energy loss of the incident parton might play an important role and 
deserve further investigation. 



1.05 



0.95 




0.85 



Figure 5: a as a function of X2, for the 158 and 400 GeV data samples. Also shown 
in the plot is a calculation of a as resulting from shadowing alone. Various shadowing 
parameterizations have been considered. 



4. Conclusions 

The NA60 experiment has measured J/ip production in p-A collisions at 
the CERN SPS, at 400 and 158 GeV, the latter energy being the lowest one 
where a detailed systematic study has been performed. The results show 
a suppression of the J/ip yield in cold nuclear matter, which is larger at 
lower incident energy. A comparison with results from previous experiments 
indicates that nuclear effects on J/tp production, at constant xf, exhibit 
a strong -y/i-dependence. The observation of different a values at the two 
energies, for a constant x 2 , m &y suggest a strong change in the J/tp break-up 
cross section and/or the presence of other effects like initial state energy loss. 
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Figure 6: Aa=a(400GcV) — a(158GcV), as a function of x 2 - The empty boxes around the 
points represent the uncorrelated systematic uncertainty. The systematic error common 
to all the points is shown as a filled box around the Aa = line. 

These data can also provide an important baseline for heavy-ion collision 
studies. 
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